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We consider two-photon exchange (TPE) in the elastic electron-proton scattering and study the 
contribution arising from the production of A(1232) resonance in the intermediate state. We calcu- 
late all three TPE amplitudes (generalized form factors), and find that the A contribution mainly 
influences generalized electric form factor (contrary to the elastic contribution, which affects mag- 
netic form factor), and the effect grows with Q 2 . If the corresponding correction is applied to the 
recent polarization transfer measurements of proton form factors, their results will change markedly. 
Thus we suggest that TPE corrections due to inelastic intermediate states are important to polar- 
CI ization experiments at high Q 2 , and should not be neglected. 
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I. INTRODUCTION 



Due to smallness of the fine structure constant a w , the elastic electron-proton scattering amplitude is dominated 
by the first-order term, corresponding to the exchange of a single photon, Fig. [TJ left. The one-photon exchange 
(OPE) amplitude has a specific structure, which allows, e.g., for the Rosenbluth separation of form factors. In the 
next (second) order the only non-trivial diagram is two-photon exchange (TPE), Fig. [TJ right. Despite its smallness, 
in some cases TPE correction is very important, because it changes qualitatively the structure of the scattering 
amplitude. Thus, TPE influence naturally explains the discrepancy between Rosenbluth and polarization methods in 
O j| proton form factor (FF) measurements. For a further review and up-to-date bibliography see e.g. Ref. 

The full TPE amplitude may be split into separate contributions according to hadronic intermediate state (IS), 
which depicted as the blob in TPE diagram. We will have the elastic contribution (pure proton IS), and inelastic ones, 
which come from multi-particle states such as pir, pmr, prj and so on. We may also distinguish the contributions of 
resonances, such as Delta, Roper etc. (This is so-called "hadronic approach". There are also QCD-style calculations 
[2L01; which assume IS to be a set of quarks. In this paper we resort to the former approach). 

At first, the TPE amplitude was approximated by the elastic contribution, which is the most well-studied one. 
Then, it was proposed to study contributions of different hadronic resonances as ISs. Kondratyuk et al. performed 
calculations with Delta Q and several other light resonances [5j . They studied the TPE correction to the cross-section 
and concluded that Delta contribution is the largest among all resonances, though still much smaller than the elastic 
one. Similar results were obtained for target normal spin asymmetry, the observable that is related to imaginary part 
y£$ ; of the TPE amplitude @. 

In this paper we present new results for Delta resonance contribution to the TPE amplitude in the elastic ep 
CN| ' scattering. To calculate it we employ the dispersion method, which was developed in Ref. Q for ep scattering and 
applied to e7r scattering in Ref. [8[. The method is described in detail in those papers, here we just recall that it 



• ensures correct behaviour of the TPE amplitude at e — > 1 (the amplitude goes to zero), 

• eliminates a need for off-shell FFs. 

In comparison to previous works of Kondratyuk et al. 0, @, we consider not just cross-section correction, but all 
three generalized FFs (TPE amplitudes) and discuss corrections to polarization transfer (PT) experiments, which was 
never done previously. 



FIG. 1: One- and two- photon exchange diagrams 
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II. THEORETICAL BACKGROUND 

To describe the elastic scattering amplitude in the presence of TPE, we will use the amplitudes Qe, Qm and Q3, 
denned in Ref. Q: 

Qe = F 1 -tF 2 +vF 3 /AM 2 

Qm = F 1 +F 2 +evF 3 /AM 2 (1) 
£ 3 = vF 3 /4M 2 

where Fi are related to the scattering amplitude as 

A-na - ( 1 P M - - \ 

M fi = —jrtw ■ u {^ F i - + m^ KF3 J u (2) 

and all notation is identical to that of Ref. 0- The TPE contribution will be indicated by the prefix S, viz. 

Qe = G e + SQ e = G e + SQ { f + SQ { E A) + ... (3) 

where Ge is usual proton FF, SQ^ is elastic contribution, SQ^ is Delta resonance contribution, and the contributions 
of other ISs (neglected hereafter), are indicated by the ellipsis. As the amplitude Q3 is absent in OPE approximation, 
it coincides with the corresponding TPE contribution, £3 = 6Q3. 

Recall that the observables are expressed via these amplitudes as follows. The reduced cross-section correction is 

Sa R = 2 Re (eG E SQ E + tG m SQ m ) (4) 

and the main contribution to it comes from SQm 0]) the correction to the FF ratio 1 R — Ge/Gm is 

SR = R Re ( - S ^ M - £ ^ Z £ ) S & 3 \ (5) 
\ Ge Gm 1 + s Gm J 

and the last term in the brackets has little effect because of small factor E ^+^ ■ 

It is interesting to note that the same combination of the amplitudes determines target normal spin asymmetry [f| : 

a Arm s 2 Q MR T ( s $e 5Q M e(i-e) sg 3 \ 

= - + 7) Q2 + AM2R2£ Im ^— - — - — j . (6) 

The AN j* vertex, in general, contains three FFs (magnetic, electric, and Coulomb ones). However, it is well- 
known, that electric and Coulomb FFs are small and experimental data are described rather well by single magnetic 
FF. Thus, to simplify the calculation, we assume purely magnetic A — » ./V7* transition, which may be described by 
the following amplitude (see e.g. fl(il|): 

Man^ = it^vmi uv a (7) 

where p and q are Delta and photon momenta, respectively (the nucleon momentum thus will be p — q) , U is nucleon 
spinor, and V a is Rarita-Schwinger wavefunction of Delta resonance, normalized according to V a V a = — 2Ma- The 
transition FF was expressed as 

i=i y l 

with m,i and c% given in Table HI These values were obtained by fitting experimental data from Ref. We will 

neglect the width of the resonance, as it was done in Ref. [![. 

For the elastic proton FFs we use parameterization from Ref. [f2[ . 



What is measured in polarization experiments is proton L/T polarization ratio, which, in OPE approximation, is proportional to Ge/Gm 
ratio. It will be convenient to divide polarization ratio by the appropriate kinematical factor and define experimentally measured FF 
ratio R, which equals to Ge/Gm m OPE approximation but really differs because of TPE corrections. 
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mi 


0. 


2.170270 


0.660810 


0.715202 


0.768494 


a 


-3.377428 


0.072839 


-20.794000 


69.497989 


-45.399399 



TABLE I: Parameters of TV — > A transition form factor, Eq. 




FIG. 2: (Color online) Delta resonance contribution to TPE amplitudes, e = 0.25, real part (a) and imaginary part (b). 



III. RESULTS AND DISCUSSION 



Figure [2] displays Q 2 dependence of Delta resonance contribution to the TPE amplitudes at fixed e (e = 0.25). As 
usual, we consider " normalized" TPE amplitudes (divided by the proton magnetic FF) . The scale is made logarithmic 
in Q 2 for better display of low-Q 2 region. 

Looking at the real parts of the amplitudes, we see sharp peaks coinciding with the resonance position. It was 
noted in Ref. Q, that such peaks are artefacts, appearing due to assumed zero resonance width. With a finite width, 
the curve must become "smeared" and the peaks should disappear. But this means that the TPE amplitude in the 
close vicinity of resonance is not adequately described by the present " zero- width" calculation. Further we will mainly 
concentrate on high-Q region (Q 2 > 1 GeV 2 ), where we do not hit the resonance and the problem will not emerge. 

We see that at high Q 2 Delta contributions grow (in absolute value) with Q 2 . Though the contribution SQj^ /Gu 
changes sign at Q 2 fa 5 GeV 2 , it still grows beyond this point (not shown in the figure). The elastic contribution 
has similar property (Fig. 0]). The difference is that the largest contribution goes to the amplitude Qe (much more 
larger than to Gm)- This fact has not much effect on the cross-section, but implies relatively large corrections to the 
polarization ratio (see below). 

The imaginary parts (Fig. [2}d), naturally, have a step-like behaviour, i.e. they vanish below the threshold and are 
non-zero above it. Having obtained the imaginary part of the amplitudes, we can perform some cross-checks of our 
results. First, we can check the sign of the TPE amplitudes with the help of optical theorem. It reads 



ImMu = 2\p\y/sa 



(9) 



where Aiu is forward scattering amplitude, s is cm. energy squared and a is total cross-section, ep — > eX. The same 
holds true for contribution of each IS h separately, i.e. 



here cr^ h > is the cross-section for ep — » eh. Putting u' — u and U' — U in Eq. ((2|), we easily get 

ImMu = 4 ^!' y Im Q E 



(10) 



(11) 



thus the optical theorem implies 



TmSQ^ 



> for Q ->• at fixed s 



(12) 
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FIG. 3: (Color online) TPE amplitudes at Q 2 = 0.5 GeV 2 (a) and 2 GeV 2 (b) 
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FIG. 5: (Color online) TPE correction to cross-section at 
FIG. 4: (Color online) Comparison of Delta (solid) and elastic Q 2 = \ GeV 2 (thin) and 3 GeV 2 (thick). Proton contribution 
(dashed) contributions, e = 0.25. (dashed), Delta contribution (dash-dotted) and total (solid). 



Note that for the elastic contribution such a check constrains only infra-red divergent part. 

We also have reproduced our results for target normal spin asymmetry from Ref. [||, using Eq. 

The e dependence of the TPE amplitudes (real parts) is shown in Fig. [3] It is substantially nonlinear for 
Q 2 < 1 GeV 2 (a) and becomes almost linear for Q 2 > 1 GeV 2 (b) for all three amplitudes. Recall that the elas- 
tic contribution SQ^ , is almost linear in e, and thanks to this fact Rosenbluth plots remain linear even with the 
corresponding correction taken into account. It is interesting to study how much nonlinearity is introduced by Delta 
contribution. We apply the method of Ref. namely, we fit calculated OPE cross-section plus TPE correction at 
fixed Q 2 by the quadratic function of e: 

<T + 5a = P [l + P 1 {e- 0.5) + P 2 (e - 0.5) 2 ] (13) 

This gives us the nonlinearity coefficient P 2 as a function of Q 2 . The Delta contribution to nonlinearity coefficient 
turns out to be rather small, |P2 A) | < 0.014 for 0.5 GeV 2 < Q 2 < 5 GeV 2 , whereas the total (elastic + Delta) 
contribution varies from 0.005 to —0.06. This should be compared to experimental value P 2 — 0.019 ± 0.027 [l3l |. 

In Fig. [S] we plot the TPE correction to the cross-section at Q 2 — 1 GeV 2 and Q 2 = 3 GeV 2 . This is the same 
quantity as in Fig. 2 of Ref. [J], and is in qualitative agreement with the latter. 

As long as now we have individual TPE amplitudes, we may easily obtain the TPE correction to FF ratio, Eq. ([5]). 
It is plotted in Fig. [6l for e — 0.5. Of course, one must keep in mind that this correction is also e dependent, as the 
TPE amplitudes are. The correction grows rapidly as Q 2 — > oo due to 
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FIG. 6: (Color online) TPE correction to measured form factor FIG 7. ( Co i or online) Results of PT experiments, with 
ratio at e = 0.5. Proton contribution (dashed), Delta contri- ( solid sym bols) and without (hollow symbols) TPE correction, 
bution (dash-dotted) and total (solid). Points are slightly offset in Q 2 for clarity. 



Expt. Q 2 e 


^Rexp ±stat . ±syst . 




3.98 0.71 
[14] 4.76 0.59 
5.56 0.45 


0.517 ±0.055 ±0.009 
0.450 ±0.052 ±0.012 
0.354 ± 0.085 ± 0.019 


0.031 
0.050 
0.078 


5.17 0.37 
[15] 6.70 0.51 
8.49 0.24 


0.443 ± 0.066 ± 0.018 
0.327 ±0.105 ±0.022 
0.138 ±0.179 ±0.043 


0.081 
0.089 
0.221 



TABLE II: Form factor ratio, measured in PT experiments (fiR e x P ) and corresponding TPE corrections (fiSR). The explicit 
factor of [i = 2.793 appears here as it is not included in our definition of R. 




• large amplitude SGe, growing with Q 2 , 

• smallness of FF ratio R itself, as it tends to zero near Q 2 ~ 10 GeV 2 . 

In Table fill the total correction, 8R = 8R^ + 8R( A \ is shown for the kinematical conditions of experiments 
[lil [To] ]. As the correction is always much larger than the quoted systematic error, it clearly needs to be taken into 
account in polarization measurements at high Q 2 . With this TPE correction applied, the FF ratio becomes negative 
already at Q 2 = 8.5 GeV 2 (Fig. 0). But is our scheme of calculation (elastic + Delta ISs) perfectly adequate for 
Q 2 ~ 5 - 10 GeV 2 ? We think that at least one needs to estimate the contributions of other prominent resonances as 
well as multi-particle states before we may apply the correction to data. Whether these contributions are small? It is 
not quite clear. Even if they are, there are many ISs that contribute, and it is not clear what will be the total effect. 

Of course, it would be nice to have a QCD-style calculation for this observable. Unfortunately, this is a hard task. 
Leading-twist QCD calculation yields only two amplitudes, Qm and Q3, since a virtual photon (gluon) cannot flip 
quark spin. The calculation of electric FF Ge or TPE amplitude 5Qe requires, at least, knowledge of quark transverse 
momenta distribution. 

Summarizing, we believe that our results give a strong indication that TPE corrections coming from the inelastic 
intermediate states may be of great importance to polarization measurements at high Q 2 , and thus deserve further 
thorough investigation. 



[1] J. Arrington, P. G. Blunden, W. Melnitchouk, Prog. Part. Nucl. Phys. 66, 782 (2011): larXiv:1105. 0951vl [nucl-th]. 

[2] Y.C. Chen, A. Afanasev, S.J. Brodsky, C.E. Carlson and M. Vanderhaeghen, Phys. Rev. Lett. 93, 122301 (2004). 

[3] D. Borisyuk, A. Kobushkin, Phys. Rev. D 79, 034001 (2009). 

[4] S. Kondratyuk, P.G. Blunden, W. Melnitchouk, J. A. Tjon, Phys. Rev. Lett. 95, 172503 (2005). 

[5] S. Kondratyuk and P.G. Blunden, Phys. Rev. C 75, 038201 (2007). 



[6] D. Borisyuk, A. Kobushkin, Phys. Rev. C 72, 035207 (2005). 

[7] D. Borisyuk, A. Kobushkin, Phys. Rev. C 78, 025208 (2008). 

[8] D. Borisyuk, A. Kobushkin, Phys. Rev. C 83, 025203 (2011). 

[9] D. Borisyuk, A. Kobushkin, Phys. Rev. C 76, 022201(R) (2007). 
[10] H. F. Jones and M. D. Scadron, Ann. Phys. 81, 1 (1973). 
[11] P. Stoler, Phys.Rev. D 44, 73 (1991). 

[12] J. Arrington, W. Melnitchouk, J. A. Tjon, Phys. Rev. C 76, 035205 (2007). 

[13] V. Tvaskis et al, Phys. Rev. C 73, 025206 (2006). 

[14] A.J.R. Puckett et al, Phys. Rev. C 85, 045203 (2012); larXiv:1102.5737H [n ucl-ex], 
[15] A.J.R. Puckett et al, Phys. Rev. Lett. 104, 242301 (2010); laxX"iv:1005.3419i rl [nucl- 



